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ABSTRACT 

First generation space-based optical coronagraphic telescopes will obtain images of 
cool gas and ice giant exoplanets around nearby stars. The albedo spectra of exoplan- 
ets lying at planet-star separations larger than about 1 AU-where an exoplanet can be 
resolved from its parent star-are dominated by reflected light to beyond 1 fim and are 
punctuated by molecular absorption features. Here we consider how exoplanet albedo 
spectra and colors vary as a function of planet-star separation, metallicity, mass, and 
observed phase for Jupiter and Neptune analogs from 0.35 to 1 /im. We model Jupiter 
analogs with lx and 3x the solar abundance of heavy elements, and Neptune analogs 
with 10 x and 30 x solar abundance of heavy elements. Our model planets orbit a solar 
analog parent star at separations of 0.8 AU, 2 AU, 5 AU, and 10 AU. We use a radiative- 
convective model to compute temperature-pressure profiles. The giant exoplanets are 
found to be cloud-free at 0.8 AU, possess H2O clouds at 2 AU, and have both NH3 
and H2O clouds at 5 AU and 10 AU. For each model planet we compute moderate 
resolution (R = A/ A A ~ 800) albedo spectra as a function of phase. We also consider 
low-resolution spectra and colors that are more consistent with early direct imaging ca- 
pabilities. As expected, the presence and vertical structure of clouds strongly influence 
the albedo spectra since cloud particles not only affect optical depth but also have highly 
directional scattering properties. Observations at different phases also probe different 
volumes of atmosphere as the source-observer geometry changes. Because the images 
of the planets themselves will be unresolved, their phase will not necessarily be imme- 
diately obvious, and multiple observations will be needed to discriminate between the 
effects of planet-star separation, metallicity, and phase on the observed albedo spectra. 
We consider the range of these combined effects on spectra and colors. For example, 
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we find that the spectral influence of clouds depends more on planet-star separation 
and hence atmospheric temperature than metallicity, and it is easier to discriminate 
between cloudy lx and 3x Jupiters than between 10 x and 30 x Neptunes. In addition 
to alkalis and methane, our Jupiter models show H2O absorption features near 0.94 
jum. While solar system giant planets are well separated by their broad-band colors, we 
find that arbitrary giant exoplanets can have a large range of possible colors and that 
color alone cannot be relied upon to characterize planet types. We also predict that 
giant exoplanets receiving greater insolation than Jupiter will exhibit higher equator to 
pole temperature gradients than are found on Jupiter and thus may exhibit differing 
atmospheric dynamics. These results are useful for future interpretation of direct imag- 
ing exoplanet observations as well as for deriving requirements and designing filters for 
optical direct imaging instrumentation. 

Subject headings: radiative transfer — scattering — methods: numerical — planets and 
satellites: general 



Introduction 



The past few years have seen the first images that resolve exoplanets from their primary stars 



most notably the t hree planets orbiting HR 8799 seen in thermal emission (jMarois et al 



2008 



Serabyn et al.ll2010l ). While more such young, hot planets are sure to be seen in the coming decade, 
space-based coronagraphs are required to achieve the planet-star contrast ratios (10~ 7 to 10~ 10 ) 
needed to detect older, cooler giant planets in reflected light in the optical (~ 0.2-1 fim). Many 
such planets will likely be observed as part of efforts to directly image terrestrial planets at high 
contrast ratios, focusing on planet-star separations close enough to be within the Habitable Zon^ 
of F, G, K stars like our G2V Sun. Here we present a range of model exoplanet albedo spectra and 
colors with the goal of improving their in terpretation and contributing to instrument design anc 



observing plans. We build on prior wo rk ( Mar lev et al 



2000 



Ford et al. 2001: Burrows et al 



1999a 



2004: 


Burrows 


2005 





Sudarskv et al 



2000 



Seager et al 



Dvudina et al.l bood : ISudarskv et a . 



200a . bood ) by developing additional model capabilities for generating albedo spectra at different 



phases and by examining a wider range of plausible atmospheric compositions, planet-star sepa- 
rations, and planet masses and radii than previously considered. In this work, we develop a tool 
compatible with our existing giant planet atmospheric radiative transfer models, and do not include 
the effects of polarization. We refer readers interest ed in the u s e of p o larimetry in the d e tectio n 
and characterizat io n of exoplanets to work such as IStam et al.1 (120041 ) ; IStam fe Hovenierl (J2005J) ; 
Stam et al.l tod ): buenzli fc Schmidl l|200g| ). 



1 Habitable Zone (HZ) generally refers to the luminosity-scaled equivalent of ~1 AU separation from our Sun; since 
our focus is on giant exoplanets, we do not require a more accurate definition of the term. 
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While numerous studies have explored the effects of system geometry and atm ospheric com 



position on exoplanet albedo spectra for trans iting close-in gas giant exoplanets (jFortnev et al 



2008bl ; iRauscher et al.ll2008l ; ISpiegel et al.ll2010l ). thermal emissions make significant contributions 
to these observati ons. Even the far red spectra of these clo s e-in exoplanets can b e dominated by 
thermal emission (jMarlev et al.l Il999al : iBurrows et al.l [2004J; iFortnev et al.l l2008al ) . In this work, 
we focus instead on understanding the reflected-light spectra of widely-separated gas and ice giant 
exoplanets that are more likely to resemble the early direct imaging targets. 

The planets themselves will be unresolved; without additional information (e.g., radial velocity 
detections) the observed phase will be uncertain. This motivates us to examine the effects of planet 
phase as well as planet-star separation, mass, and metallicity on the reflected-light albedo spectra 
of Jupiter and Neptune analogs. We use the term 'metallicity' to refer to an enhancement in heavy 
elements compared to solar elemental abundances. 



1.1. Contrast and phase 

The contrast ratio between an exoplanet and its star depends on many factors, including 
the physical properties of the planet and star, the system geometry, background sources, and the 
instrumentation used. A simple way of estimating the required contrast is to assume that the 
planet reflects its incident starlight isotropically, such that its apparent brightness is const ant over 
its ill uminated surface, as a Lambertian surface. Then the contrast can be expressed as (jSobolev 



19751 ): 



C („) = *A. W (£) 

where a is the planet phase angle, R p is the planet radius in km, d is the planet-star separation 
in km, and A g is the geometric albedo, which we formally define in §2.21 The geometric albedo 
generally takes values between and 1 for the fraction of monochromatic light the planet reflects 
towards the observer at full phase, although it can be larger than 1 for anisotropic scattering 
atmospheres or surfaces. If one assumes the planet is at quadrature, a = vr/2, and that the 
geometric albedo at the wavelength of interest is 0.5, then for a Jupiter analog (lRj), the contrast 
at 0.8 AU is ~ 3.8 x 10~ 8 and at 10 AU is ~ 2.4 x 10 -10 . For a Neptune analog (1 R N ), the contrast 
at 0.8 AU is ~ 4.5 x 10~ 9 and at 10 AU is ~ 3 x HT 11 . For an Earth analog (1 Re) at 1 AU with 
A g = 0.3, the contrast is ~ 1.2 x 10~ 10 . High contrast is more difficult to achieve closer to the 
star than further away, however, the d~ 2 dependence also means that more distant objects will be 
fainter in reflected light. 



sin a + (tt — a) cos a 



IT 



(1) 
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1.2. Composition 



As discussed in lMarlev et al.l (|2007l ). determining atmospheric composition is a central objec- 
tive for characterizing extrasolar giant planets. A census of how composition varies with planet 
mass and orbital distance will help inform theories of giant planet formation. The giant planets 
in our solar syste m motivate us to consi der a large range of metallicities in our model exoplanets. 
As summarized in lFortney et al.1 ( 2008a ). the Galileo and Cassini spacecraft have shown that the 
atmospheres of Jupiter and Saturn have enhanced heavy element abundances compared to that 
of the Sun. For example, Jupiter has ~ 2x to 4x s olar abundances of oxygen, carbon, nitro- 
gen, sulfur, and various noble gases (lAtreya et al.ll2003l ). Understanding how giant planets formed 
and the role of processes such as plane te simal bombardment and accumulation durin g formation 



(jOwen et al.l 1 19991 ; iGautier et al 



2001a 



bllGuiilot & GladmarJ[200ol : klibert et alJbood ). erosion of 



the heavy-element core (|Stevensonlll985l : lGuillot et al.ll2004l ). direct ac cretion of metal-rich disk gas 
(IGuillot fc Huesol 120061 ) . or chemical fractionation within the planet (jStevenson &: Sarpeterlll977l ; 
Lodderal2004l ) in enriching gas giants in heavy elements is relevant both to understanding our own 



solar system as well as exoplanet systems. 

This motivates our goal of better understanding the combined effect of planet-star separation, 
composition and planet phase on spectra that are likely to be obtainable in early reflected-light 
images of exoplanets. For this initial study, we assume simple circular orbits where the planet-star 
separation is constant, and examine how the albedo spectra and colors change with planet phase 
due to the change in angle from the source to the observer. More complex orbits, such as eccentric 
orbits, require consideration of how the structure of the exoplanet's atmosphere would evolve as a 
function of distance from its parent star. An initial study of the effect of eccentric orbits, using 
interpolated temperature and pressure profiles correspond ing to different planet -star separations, 
was presented for a lMj extrasolar giant planet (EGP) in Sudarsky et al. (|2005l ). 



1.3. Clouds 



As with any planet, exoplanet albedo spectra are very sensitive to the presence and structure of 



clouds (jMarlev et al.lll999al : 



Sudarsky et al.ll2000i ). Becaus e the global he i ght an d optical thickness 



of cloud layers are difficult to model on an a priori basis, iMarlev et al.l (| 1999a! ) c onsidered a few 



charac teristic cloud models in order to test the sensitivity of albedos to clouds. ISudarskv et al 



(|2000i ) likewise employed a fairly simp le method for computing c loud opacity. In this study, we 
use the more complex cloud model of lAckerman &: Mar ley! (|200ll ) to calculate the height, parti- 
cle sizes and consequently optical thickness of water and ammonia clouds. A comparison of the 



m 



Ackcrm an fc: Mar lev! (1200 if) m odel with the models used inlMarlev et al.l (Il999al ) can also be found 



Ackerman &: Marleyl (|200ll ). The lAckerman &: Marleyl (|200lf ) cloud model provides wavelength 



dependent measures of the anisotropic scattering in terms of the asymmetry factor g, and single- 
scattering albedo Co of the cloud particles at each temperature and pressure level. In 53 we discuss 
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how these parameters are used to approximate anisotropic scattering as well as describe planned 
modifications to the cloud model to provide additional scattering information. 



The ISudarskv et al.l (|2000l . 120051 ) models generally focused on 1 Mj planets and considered 
the effects that might be manifested in planets with inclined or highly elliptical orbits (such as 
transitioning from cloudy to cloud- free and thus bright to dark in reflected light). The effect of 
geometry on the orbit determ i nation of close-in exoplanets from p hotometry for was als o recently 
presented in detail by iBrownl (|2009l ). Here we exte nd the work of Marlev et al.l (|1999al ) to model 



exoplanet observations as a function of phase and ISudarskv et al 



Lev 

(|200a l2005h by considering a 



greater range of planet masses and metallicities. In this work we consider both cloud-free (0.8 AU) 
and cloudy (> 1 AU) planets. 



1.4. Organization 

We briefly review the basics of exoplanet albedo theory in §[2] as framework for the following 
results and discussion. While this material is familiar to the solar system planetary science com- 
munity, it bears reinforcing because of the subtleties between various types of planetary albedos 
and their relationship to phase. In $3] we describe how we generate a series of exoplanet model at- 



model that was previously tailored for use with EGPs ( 


Marlev 1997; Fortnev et al. 


2005. 


2006; 


Marlev et al.ll2007 


Fortnev Sz Marlev 


2007 


Fortnev et al.ll2008a 


b ) and brown dwarfs ( 


Marlev et al. 


1996; 


Burrows et al. 


1997; 


Marlev 


1997; 


Marlev et al. 


2002; 


Saumon et al. 2006, 


2007) 


. The model 



was originally used for solar system planetary bod ies such as Titan and Uranus (jToon et al.lll977l . 
19891 : McKay et al.lll989l : Marlev k McKay! Il999bl ). These global mean ID exoplanet models are 
then used as input to a high resolution albedo spectra model. We describe our updates to the 
albedo model that allow us to calculate emergent intensities and thus albedo spectra as a function 
of phase. 

In 2] we present results: albedo spectra from 0.35 /iin to 1 /jm as a function of planet type, 
planet-star separation, metallicity, and phase. Because of the practical constraints inherent in 
spacecraft flybys of solar system giant planets, there is relatively little data in the literature report- 
ing planet- averaged phase functions; we do compare Voyager 1 data of Uranus and Neptune from 
Pollack et al.1 (| 19861 ) with our model phase func tions. We a l so com pare the albedo spectra with 
observations of solar system giant planets from [Karkoschkal ( 1994 ). Since early exoplanet direct 
imaging observations will have limited resolution, in £j5]we consider lower-resolution, coarse spectra 
derived from our high-resolution results. We consider R = A/AA = 5 and 15, and we also consider 
how changes in albedo spectra might present themselves in terms of color-color comparisons using 
standard filters in the optical. This helps us to evaluate features that might be detectable using a 
combination of different wide (or narrow) filter bands. 



Color-color diagrams were suggested as a comparative analytic tool for direct imaging of exo- 
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planets i n an example that con sidered the color diversity of solar system planets in iTraubl (|2003l ). 
Recently, iFortnev et al.l (I2008al) presented d etailed color-color comparisons of hot young Jupiters in 
the infrared. Earlier, ISudarskv et al.l (|2005l ) modeled reflected light albedo spectra and performed 
color-color analyses in the optical of a 1 Mj EGP using a similar approach but with a different 
model and implementation than that used in this work. While our 1 Mj results tha t include clouds 



tend to be less red at full phase and less blue at new compared with those in ISudarskv et al 



(|2005l ). we generally agree with their cloud- free result and extend their approach to include dif- 
ferent compositions, Neptune analogs, and comparisons with observed spectra and colors of solar 
system planets. 



2. Background 

This section provides some background on exoplanet direct imaging methods and the scientific 
value of albedo spectra constructed from observations made in the optical. We briefly review 
common direct imaging terminology and techniques in E j2.ll and then formally define geometric 
albedo and illustrate how albedo spectra of exoplanets relate to exoplanet science goals in §2.21 
We also briefly discuss the challenges and current limitations that affect direct imaging efforts. 
Consideration of both the instrumentation constraints for direct imaging as well as exoplanet science 
goals help to define and justify the planet-star separations and atmospheric model types used in 
this work. 



2.1. Direct Imaging 



To give the reader an idea of the challenges inherent in coronagraphic detection and characteri- 
zation of gas or ice giant exoplanets, Figure [1] shows simulated direct images of Jupiter and Neptune 
analogs observed in t hree different 100 nm wide bands around a solar analog at a distance of 10 
pc from the observer (jCahoy et al.ll2009l ). The simulations are performed for a 4 m diameter space 
telescope, where the light from the parent star has been suppressed by a Phase Induced Amplitude 
Apodization (PIAA) coronagraph (IGuyon et al.ll2005l ). and stellar leakage is included as one of 
several background sources. Importantly, for this figure only, we assume a gray albedo spectrum 
with a geometric albedo that is a constant 0.3 with wavelength and thus the same for e ach planet 
in each bandpass. The simulations are described in greater detail in lCahov et al.l (|2009l ). 



Figure [T] helps to emphasize the importance of planet-star contrast ratios and angular res- 
olution in a direct imaging observation. It also puts into context several factors that are not 
specific to a particular instrument or method, but that are important to consider when both in- 
terpreting observations and when constructing exoplanet models that correspond to early direct 
imaging capabilities. For example, the simulated images make clear the relevance of planet size, 
wavelength-dependent brightness (albedo spectra, one can imagine the effect of a decreasing albedo 
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with wavelength on the images), minimum and maximum detectable separations from the parent 
star, contributions from background sources like exozodiacal debris disks, and the effect of system 
geometry and planet phase on the observations. The wavelength dependent plate scales in this ex- 
ample also give a sense of the level of calibration required to construct coarse albedo spectra even 
from concurrent measurements. For photon-starved exoplanet observations, the tradeoff between 
resolution, integration times, and noise is the major challenge in obtaining data in multiple bands. 

The term inner working angle (IWA) is often used to describe the minimum detectable sepa- 
ration between a star and planet in a high-contrast direct imaging measurement, where the IWA is 
generally substantially larger than the diffraction limit (e.g., the 1.22A/-D Rayleigh criterion). For a 
band-limited coronagraph, the IWA is defined as the separation at which the planet's light is at half 
of the ma ximal throughput of the starlight suppression technique at a contrast level sufficient for 
detection (ICagigal et al.ll2009l ). An instrument that is improved to achieve a smaller IWA can de- 
tect exoplanets at closer star-planet separations and thus can also detect exoplanets around stars at 
greater distances from the observ er. The IWA for curren t observing systems is typically on the order 
of a few tenths of an arcsecond (jBeichman et al.ll2010l ): instruments in development may achieve 
IWAs on the order of a tenth of an arcs econd. Examples of som e current direct imaging methods 



include high-performance cqronagraphs (IKuchner Traub 



2005 



Trauger fc Traubll2007l: 



(Wallace et al 



used by 



200 



0; 



Mawet et a 



Martin Levine et al 



Marois et al 



2002 



Kasdin et al 



2003 



Guvon et al 



(2006) and lLafreniere et all (j2007h . 



20Qgh. external occulters (jCashlbood l. and optical nullers 
20o3L 

as well as post-processing algorithms such as those 



In practice, for current direct imaging methods the IWA is typically larger than 2X/D at the 
required contrast levels for imaging exoplanets in reflected light (the IWA for the PIAA coronagraph 
simulated in Fig.Q]is also ~ 2X/D). Early optical direct images of exoplanets will likely be limited to 
planet-star separations greater than 2X/D, since the IWA metric generally does not include external 
contributing factors such as noise, background sources, coatings, optical throughput, quantum 
efficiency, etc. There are also platform-specific limitations, for example, primary mirror diameter 
D is limited on space-based direct imaging missions by either deployment complexity, cost, or launch 
vehicle fairing size. Ground-based instruments do not have this limitation on D, however, they do 
have other challenges, such as contrast limitations due to turbulence in Earth's atmosphere. Given 
the capability of current technology, using a minimum separation of ~ 1 AU for our exoplanet 
models around a Sun-like star seems reasonable for this first set of simulations. In addition, a 
separation of ~ 1 AU is far enough away from a Sun -like star that reflect ed light will dominate the 
features in the optical spectrum, even of gas giants (jMarlev et al.lll999al ). 



The model planets shown in Fig. [T] are based on their solar system counterparts and are 
separated from their parent star by luminosity-scaled equivalents of 1 AU, 2 AU, and 5 AU. We 
hold the albedos constant vs. A at 0.3 in these simulations, to better illustrate the other wavelength- 
dependent changes. If we were to ignore the wavelength-dependent changes in albedo, the larger- 
radius gas giant planets are simply brighter and easier to detect than terrestrial exoplanets. This 
is particularly true at redder wavelengths where the A dependence of the IWA 'pushes' the planet 
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closer to the IWA, effectively making smaller separations harder to detect at longer wavelengths. 
As we discuss in £JH realistic giant planets get darker at smaller planet-star separations when the 
temperatures become too warm for bright clouds to form. Depending on the star, the darkening 
at temperatures that are too high for cloud formation may occur inside of the IWA, making the 
detection of such planets less likely. For reflected- light direct imaging of cooler exoplanets, Fig. Q] 
show that it is easier to detect the larger Jupiter and Neptune analogs compared with terrestrial 
planets, which motivates our use of them as the exoplanet models in this work. 

For some direct imaging methods, there may also be limitations on the maximum planet-star 
separation, or outer working angle (OWA). This can be due to the optical design or elements such 
as deformable mirrors, detector size, sub-apertures, etc. It is desirable to have an OWA of ~ 20X/D 
or larger, if possible, so that in addition to exoplanets, exozodical dust and circumstellar debris 
disk structure can be characterized at much wider separations. 

The existence of giant planets at separations large r than ^ 10 AU i s difficult to account for 



in st andard core accretion models ([Pollack et al.l Il996l ; llda &: Linl 120051 ; iDodson-Robinson et al 



2009). Their presence at very large planet-star separation s may be due to a different fo rmation 
mechanism, such as gravitational fragmentation in the disk (|Bossll2000l : iKratter et al.ll2010i ) or due 
to a combination of multiple effects that move planets formed in dense inne r regions of a disk out 
to distant orbits, including outward migration and planet-planet scattering (jVeras et al.1 12009). 



In this work, we model Jupiters and Neptunes out to 10 AU based on the theoretical predictions 
for giant planet separations and the fact that the contrast requirement for cool giant planets beyond 
10 AU is challenging (d~ 2 dependence on the incident light, as noted in fll.ip . Using Eq. Q3 a lRj 
planet with geometric albedo of 0.5 has a contrast of ~ 1 x 10~ 10 at 15 AU and ~ 6 x 10~ n at 
20 AU when observed in reflected light at quadrature. Such contrasts are comparable to or fainter 
than that expected for terrestrial planets at 1 AU, thus motivating our outer model boundary at 
10 AU. In $3] we discuss the role that planet-star separation plays in controlling the atmospheric 
structure (clouds) and albedo spectra of our model exoplanets. 

Simultaneous images taken with different filters can function as coarse spectra. In Fig. [2]we use 
observed spectra of the outer solar system planets and Titan from iKarkoschkal ([1994) to illustrate 
how coarse spectra from an early direct imaging observation might be used to differentiate between 
planet types. The original spectral data from IKarkoschkal (|1994l ) are binned down to resolution 
R = A/AA = 5 and R = 15. Even in the R = 5 spectra, one can distinguish between the ice giants 
(Neptune, Uranus) and the gas giants (Jupiter, Saturn). The ice giants are brightest in the blue 
and get darker in the red, and the gas giants are not as bright in the blue, due partly to the effects 
of photochemical products, but are brighter mid-band and in the red. An important goal of our 
study is to understand the extent to which such distinctions hold in general for exoplanets and how 
the phase at which they are observed affects the colors. 

As we discuss in more detail in §2.21 and $11 coarse spectra can be used in combination with 
atmospheric models to constrain the temperature, composition, the presence or absence of clouds 
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in an e xoplanet's atmosphere, and the exoplan et's radius (jMarlev et al.lll999al ; iMarlev fe McKay 



1999b : 



Marlev et al. 



2007: 



Fortnev et al 



2008a 



Color-color comparisons of exoplanets can be 
made (jTraubl 120031 : ISudarskv et al.l 120051 : iFortnev et al.l l2008al ). and understanding how system 



geometry and explanet composition relate to variability in the observations is important for inter- 
preting these color-color comparisons. 



2.2. Albedos 

The definition and relationships between between geometric albedo A g , phase a, the phase 
integral q, spherical albedo A s , and Bond albedo Ab, and additional parameters of interest such as 
planetary effective temperature are important to understand in the context of interpreting direct 
imaging observations. Although this material is not new, we review it since th e distinctions between 
the various types of albedos are important to emphasize. Note that while IMarlev et al.l (|1999a) 
presented geometric albedo spectra and tables of Bond albedos for exoplanets in systems with 
varying stellar types, our focus here is more closely on geometric albedo, albedo spectra, and 
planet phase functions. We do tabulate Bond albedos for our exoplanet models, but only for a 
solar analog parent star. 

Measurements of the planet's brightness at a given wavelength and full phase yields the planet's 
monochromatic geometric albedo: 



A g (A) 



F p (A,a = 0°) 
F ,l (A) 



(2) 



Geometric albedo A g (A) is the ratio of the reflected flux F p (A, a = 0°) of an object at full phase 
(a = 0°) to the flux from a perfect Lambert disk, i^^of the same radius R p under the same incident 
flux Fq at the same distance d from the star. A planet's phase angle a is the angle between the 
incident ray from the star to the planet and the reflected ray from the planet to the observer. The 
scattering angle is related to the phase angle via a = ir — G. 

In this work, we use the term geometric albedo to refer to albedo spectra at full phase. We use 
the term albedo spectra to refer to spectra observed at different orbital phases. 

In order to reduce confusion in later discussions of planet phase functions and particle scattering 
functions, we use the notation A g (A) for geome tric albe d o her e instead of using the other common 
notation for geometric albedo, p. As noted in ISobolevI (119751 ). for a perfectly reflecting Lambert 
sphere, t he geometric albedo is 2/3, an d for a semi-infinite purely Rayleigh scattering atmosphere, 
it is 3/4 (|Dlugach fc Yanovitskiilll974h . 



Flux measurements at different planet phase angles can be used to define the phase function 
$(A,a): 
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Fp (A, a) 
F (A) 



A g (X)(^j\(X,a) (3) 



As noted for Eq. [21 A g (A) is defined at a = 0°, F p (A, a) is the monochromatic planet flux, F (A) 
the monochromatic stellar flux, R p the planet's radius, d the planet-star separation, and $(A,a) 
the planet's phase function at angle a. The phase function is normalized to be 1.0 at full phase. 
The scattering phase function of the particles in the exoplanet's atmosphere, p (0) contribute to 
the disk-integrated phase function. The observed phase function allows calculation of the phase 
integral, q: 

f'Tl 

q(X) = 2 / $(A, a) sin ada (4) 
Jo 

The spherical albedo, 



A s (X) = q(X)A g (A) 



(5) 



is the fraction of incident light reflected towards all angles. 

The planet's Bond albedo, an incident flux-weighted and wavelength-integrated function of 
A s (A) and $(A, a), normalized by the incident stellar flux, can then be computed: 



A B 



f o °°A s (A)F (A)dA 
/ o °°F (A)dA 



(6) 



It is important to appreciate that since the Bond albedo is weighted by the incident flux, two 
identical planets with identical geometric albedo spectra will have di fferent Bond albedos if the 
incident flux is different (e.g., comparing an F star and an M star) (jMarlev et al.l Il999al ). The 
Bond albedo can be used to find the effective temperature, T e g of a rapidly rotating planet: 



4nRlaTj s = E int + ttR 2 p [f (1 - A B )] ^ 



(7) 



where a is the Steffan-Boltzmann constant, ttS is the incident flux of the parent star at a planet-star 
separation of 1 AU, d is the pla net-star separation, h ere in AU, and E{ n t represents the internal 
energy sources of the planet (e.g. IConrath et al.lll989l ). 



An additional redistribution parameter / may also be introduced to parameterize the global 
redistribution of energy to solve for a day and nightside temperature. For solar system giant planets, 
typically / = 1, since Jupiter or Neptune-like planets with rotation periods of hours and effective 
temperatures below 500 K have sufficiently long radiative time constants such that atmospheric 
dynamics likely to redistribute incident flux around the planet. For exoplanets, sometimes different 
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values of / are used due to the diversity of exoplanet system geometries. For example, f = 2 
represents the case where the redistribution of the absorbed flux is over the day side only, yielding 
a higher day side effective temperature. Here, we use / = 1. 

In <J3] and Appendix [A] we describe updates to our albedo spectra model that allow us to 
generate albedo spectra at different planet phases for a range of input model gas and ice giant 
exoplanet atmospheres, which can then be used to calculate exoplanet colors, phase integrals, 
spherical albedos, and Bond albedos. 



3. Description of Methods and Models 



In this section, we describe the two-part approach to generating our results. For each model 
planet we first compute a mean thermal structure that accounts for the deposition of incident 
flux and an assumed internal heatflow. The ID mean model is computed utilizing our radiative- 
convective equilibrium model accounting for transport of thermal emitted flux and incident light. 
Once we have the mean global structure in hand we compute the detailed scattering calculations 
as a function of orbital phase using our albedo model. 

In this section, we describe the two-part approach to generating our results. For each model 
planet we rst compute a mean thermal structure that accounts for the deposition of incident ux and 
an assumed internal heatow. The ID mean model is computed utilizing our radiative-convective 
equilibrium model accounting for transport of thermal emitted ux and incident light. Because 
radiative timescales are longer than a r otation period and atmo spheric and internal dynamics tends 
to distribute energy around the planet (jlngersoll fc Porcdl 19781 ). the ID mean profile is expected to 
be an excellent description of the atmospheric profile at most points on the planet q Once we have 
the mean global structure in hand we compute the detailed scattering calculations as a function 
of orbital phase using our albedo model. The ID model alone could not be used for this latter 
application since it treats mean conditions encountered at a xed stellar zenith angle. Our two step 
approach allows detailed modeling of the scattered light at arbitrary phase and spectral resolution 
and follows our sim i lar approaches at Titan a nd Uranus with earlier versions of the same models 
(jMcKay et al.lll989l : iMarley McKay! Il999bl ). For clarity we will use the description 'radiative- 
convective model' to refer to the model employed to generate the ID exoplanet models and the 
description 'albedo model' to refer to the model used to generate the albedo spectra. 



2 Note the discussion in H4.4. 1 1 speculating that some giant planets may have cold polar regions; such a temperature 
pattern would break our assumption of a homogeneous global thermal profile. 
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3.1. Radiative-Convective Exoplanet Model 



We generate the exoplanet model atmosph e res with the 1-D iterat i ve radiative-convective atmo - 
sphere model previously u s ed inlMarlevi (119971'); iBurrows et all (119971'); iMarlev et al.l (11999a! . 120021 ) 



diative tra nsfer method was develope d by iToon et al.1 (119771 . Il989l ) , imple mented bvlMcKav et al 



Fortnev et aD ibood . bood . boOTT) : iFortnev fc Mar lev! (boOTn and ll-oitnev et all (bo08allbK The 



ra- 



ti 19891 ) and Marlev & McKavi (11999b! ') and used in application to Titan (McKay et al.lll989l ) and 
Uranus (jlVIarlev fc McKavlll999bn ~ as well as to EGPs and brown dwarfs, as noted in £JU The model 
computes the deposition of incident flux from 0.268 fim to 6.064 /jm and accounts for thermal emis- 
sion out to 227 /im. 



In this work, we use elemental abundances fromlLodderd (120031') . We employ tabl e s of c h emica l 
equilibrium compositio ns following iFeglev h Loddersl (|1994i ) and lLodders fc Feglevl ([2002], 120061 ) ; 
Freedman et al.l (|2008l ). The chemistry calculations include 'rainout,' where re fractory specie s 
are depleted from the atmosphere due to their condensation into cloud decks (jLoddersI Il999l ) . 
The spectra of both brown dwarfs and our solar syst em's giant planets can only be reproduced 
when chemistry calc ulations incorporate this process (IFeglev &: Loddersl Il994l : IMarlev et al.ll2002l : 



Burrows et al. 



20021 '). 



As noted in £ 11.11 the four gas giants in our solar system show an enhancement of heavy elements 
compared to solar abundances. The amount of enhancement provides additional clues about planet 
formation and evolution. Our models include a range of compositions, with metallicities of 1 x and 
3x solar for the model Jupiters, and 10 x and 30 x solar for the model Neptunes. To be specific, 
we use metallicities of [M/H] /[M/H cosm i c ] = 0.0, 0.5, 1.0, and 1.5. So although we refer to them as 
3x and 30x, the 3x and 30x metallicities are actually 3.16x and 31. 6x. 

The radiative-convectiv e model uses the correlated- & method for the tabulation of gaseous 
opacit ies (jGoody et al The extensive opacity database we use is fully described in lFreedman et al 

(|2008l ). We note that we have yet to inc orporate the recent tabulation of methane opacities pre- 
sented in iKarkoschka fc Tomasko I (|2010j). We also note the general difficulties in determining the 
accuracy of opacities at the temperature and pressures of interest, and of reconciling laboratory 
observations with Earth-based, space-based, and in-situ observations. 



We use the cloud model of lAckerman &: Marlevi (|200ll ) to describe the location, vertical dis- 
tribution, and particle sizes of the major cloud-forming species, which are H2O and NH3 in this 
work. Clouds are coupled iteratively into the radiative-convective model using the approach of 
Ackerman Marleyl (120011 ) as described in §1.31 The model computes the vertical variation of 
particle sizes and number densities given an assumed cloud sedimentation efficiency, / se d which 
we discuss further below. Since the cloud model is converged along with the atmospheric thermal 
structure, the clouds are fully self consistent with the atmospheric thermal structure of each model 
planet. Note that we neglect NH4SH, which m ay produce an in termediate cloud between the water 
and ammonia clouds in Jupiter's atmosphere ([West et al.ll2004i ). Even at Jupiter the r elative roles 



that NH3 and NH4SH clouds play in influencing the albedo spectra remain uncertain (jlrwin et al 
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20051 ) and thus we treat only the former cloud. 



For each model planet our ID modeling thus produces a vertical temperature and pressure, 
T(P), profile, a cloud structure model, and the variation in composition of all major species. We 
use these profiles and cloud information (wavelength-dependent asymmetry factor g and single- 
scattering albedo ui at each temperature and pressure level) as input to the albedo spectral model 
described in £)3,2i 



3.1.1. Model Temperature and Pressure Profiles 

Table Q] summarizes the exoplanet model gas giants (Jupiters) and ice giants (Neptunes) used 
in this work. The internal heat flows assumed for the Jupiters and Neptunes correspond to effective 
temperatures (Ti nt ) of 100 K and 50 K, respectively, for an isolated object, and an age of ~4.5 Gyr. 
Gravity for Jupiters is 25 ms~ 2 and for Neptunes is 10 ms" 2 . We compute radiative-convective 
equilibrium T(P) profiles and low-resolution spectra for exoplanet models at planet-star separations 
of 0.8 AU, 2 AU, 5 AU, and 10 AU at lx and 3x solar metallicity for the Jupiters and 10 x and 
30 x solar metallicity for the Neptunes. As noted in §2.1} the maximum and minimum planet-star 
separations were chosen to span the range of exoplanets detectable with early coronagraphic direct 
imaging methods. 

Figure [3] shows the resulting T(P) profiles for the models. Solid lines are used for the lower 
metallicity lx and 10 x cases, and dashed lines are used for the higher metallicity 3x and 30 x 
cases. Condensation curves, courtesy K. Lodders, are shown on each plot for both H2O (warmer) 
and NH3 (cooler). The condensation curves correspond to lx (solid line) and 3x (dashed line) 
solar metallicity for the Jupiter models and 10 x (solid line) and 30 x (dashed line) solar metallicity 
for the Neptune models. 

To illustrate the impact of clouds on the albedo spectra and thus the detectability of exoplanets, 
we chose planet models across the full range of planet-star separations likely to be probed by early 
direct imaging observations. Our planet-star separations include cloud-free, H2O, and combined 
H2O and NH3 cloud cases for both the Jupiters and Neptunes. However, for the purposes of 
comparison and simplicity, we also used the same planet-star separations for both the Jupiters and 
the Neptunes to limit the number of models to a few representative cases. 



While the lAckerman fe Marleyl (120011 ) approach results in clouds with varying particle sizes as 



a function of height through the atmo sphere, for heuristic purpo ses the mean or 'effective' cloud 



particle radii can also be computed (see lAckerman &: Marleyl (|200lh ). Table[2]lists the effective radii 



of the condensate particles for the clouds in each of the exoplanet models. We stress that these 
are provided here for comparison between the models; the actual radiative transfer calculations 
do not employ these mean particle sizes, but rather the actual computed particle size profiles. 
We used the same values for f se d, the ratio of the microphysical sedimentation flux to the eddy 
sedimentation flux, for both Jupiter and Neptune models at a given planet-star separation. At 2 
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AU, we used / se d = 6, and at 5 AU and 10 AU, we used / se( j = 10. These values of / se( j were 
selected partly from our experience with fitting L dwarf spectra ((Stephens et al. M2009I ) and partly 
to allow convergence (very thick clouds modeled with low values of / se d do not converge well). 
Ackerman Marlevl (|200ll ) found good agreement with the NH3 cloud in Jupiter's atmosphere 
using = 3. Larger values of / sp h lead to clouds of slightly smaller vertical extent and optical 



depth (j Ackerman &; Marlevl l200ll ). In a future study, we will explore the effect of changing cloud 
thickness on planet colors. For the current study, all of the clouds have large optical depths, so 
the impact of varying / se( j is small. We stress that for this study we are primarily interested in the 
macro effect of clouds and not in detailed changes to the cloud model. A higher value of / se d was 
required for the water cloud planets in order to produce properly converged models. 

Fig.[3]shows that at 0.8 AU, both the lx and 3x Jupiters are just warm enough to be cloud-free. 
At 2 AU the Jupiters are just warm enough to only condense H2O clouds, and at 5 AU and 10 AU 
the Jupiters have both H2O and NH3 clouds. Fig.[3]also shows that the H2O and NH3 condensation 
curves for the 10 x and 30 x Neptunes are a bit warmer than their lx and 3x counterparts. The 
0.8 AU Neptunes are also just warm enough to be cloud-free. The 2 AU Neptunes may be just 
cool enough to condense both H2O and NH3. However, in this work we neglect the condensation of 
NH3 at 2 AU, for better comparison with the 2 AU Jupiter models that only condense H2O. The 
5 AU and 10 AU Neptune models have both H2O and NH3 clouds. 



Our lx Jupiter models appear to be a bit cooler than those shown in ISudarsky et al.1 (|2005l ). 
This is probably due to our computing a global mean temperature profile which assumes that 
the absorbed inciden t flux is redistributed by atmospheric dynamics evenly around the planet. 
Sudarskv et al.l (|2005l ) instead compute a mean dayside temperature profile without energy redis- 
tribution. While the efficiency of energy redistribution is an interesting topic for the very hot 
Jupiters, at the lower effective temperatures relevant to this study the radiative time constants 
are long and the incident flux is likely to be evenly distributed, as it is on the solar system jovian 
planets, which do not exhibit any significant diurnal thermal differences. 



3.2. Albedo Spectra Model 



Given the atmosphere models generated by the radiative-convective model we next apply the 
albedo spectra model. This albedo spectra mod e l has previously been used f or a y ariet y of planetary 



objects , and is described in detail iruToon et al.l (119771 . ll989l ): lMcKay et al.l (| 19891 ) and lMarley et al 



(|1999al ): iMarlev &; McKavl (|1999bl ). In this section we describe the configuration of the albedo 



spectra model used here. The albedo model takes as input the exoplanet model's gravity and 
atmospheric temperature, pressure, composition, and cloud properties that are generated within 
the radiative-convective model discussed in ^3.11 



The albedo spectra model uses the same chemistry and opacity database as were used for 
obtaining the thermal structure. In the thermal structure model, opacities were binned into intervals 
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using k-coefficients. For the albedo calculation we use monochromatic opacities over the ~ 790 
points from 0.35 /im to 1 //m. The opacity of course varies with altitude through the model 
and we account for the changing gaseous composition as well as the influence of temperature and 
pressure on the opacities themselves. In this work, we do not mode l phot ochemistry or attempt to 
account for photochemical hazes or chromophores (see I West et al.l (|2004l ) for a review) , which are 
responsible for the low albedos for the real Jupiter in the blue and UV (< 0.4 jjcm) and also for its 
ruddy complexion; this is discussed in more detail in iMarley et al.l (|1999al ). While we expect that 
photochemistry will play an important role in exoplanet spectra, we reserve such considerations for 
future work. 



3.2.1. Radiative transfer description 

Figure 2] shows a schematic of the plane-parallel model geometry with incident angle /Uo = 
cos/3 = sin 7] cos (£ — a) and observed angle [i\ = cos$ = sin r\ cos £, where £ and r) are planetary 
coordinates of longitude and co-latitude, respectively. 

We use a finite plane-parallel atmosphere divided into N levels and N-l layers. The plane- 
parallel framework is shown in Figure 01 Emergent intensity is I (t, fii) with optical depth n from 
level i to the top of the atmosphere and incremental change in optical depth of 5 n = Tj+i — Tj 
between levels. Both incident and emergent angles are positive from the surface normal, with 
incident light from direction fiQ = cos/3 and the observer direction [i\ = cosfi. A negative sign is 
used to denote the direction of incident flux, e.g., — hq. With a source at fio, the exact expression 
for the emergent intensities through the N levels is: 

/(T i , W ,0i) = J(r i+ i,//i,0i)e- ,5 -> 1 + [ Ti Stfinx^e-'fn— (8) 

Jo IJ-i 

where S (r',yUi, (f>i) is the source function and contains contributions from both direct and diffuse 
scattering. The incident collimated flux crosses the unit area perpendicular to direction defined by 
— an d u) is the 'single-scattering' albedo, a parameter used to represent the amount of scattering 
taking place versus absorption. u> is calculated for each layer and at each wavelength by taking the 
ratio of the scattering coefficients to the total scattering and absorption coefficients in a layer. The 
asymmetry parameter, g, is similarly calculated for each layer and wavelength using asymmetry 
contributions and opacities of the layer constituents. The scattering and asymmetry contributions 
include an approximation to Raman scattering, Rayleigh scattering and clouds in this model. 

While we begin with an exact solution to the radiative transfer in Eq. (8), we introduce an 
approximation by substituting a two-stream solution for the scattered radiation field (intensity) in 
the source function. In the case of negligible diffuse radiation (when single scattering dominates, as 
it tends to for clear atmospheres) this approximation has no effect. As the importance of multiple 
scattering increases (e.g., when there are clouds and hazes), the influence of this approximation 
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grows. We present our two stream solution with so me detail, but str ess that the albedo calculation 
as used here is not simply a 'two stream solution' (jToon et al.lll989l ). 



The radiative transfer equation is often solved using approaches that do not preserve the 
a ngular dependence of the s cattering phase function, su ch as th e two-s tream methods summarized 



m 



Meador fc Weaver! (ll98Gh or th e Feau trier method in iMihilasI (|1978l ). Here we follow the 'source 



function' approach in iToon et al.l (|1989l ) where we first use two-stream quadrature to solve for the 
diffuse scattered radiation field. We use the resulting two-st ream intens i ty It a s an approximation 
to the source function (S v versus S v t in Eq. 1 and Eq. 52 of IToon et al.l (|1989i )). We are then able 
to approximate the angular dependence (p,<j>) and (//,<?/) in the scattering phase function. We 
use the two-stream 'source function' approach here largely because of its compatibility with the 
existing model. In a future paper we will compare the approach in this work with other radiative 
transfer methods that pres erve the angular dependence of the scattering phase function, such as 



the 5— M stream method (jWiscombe 1 119770. Compa ri sons with models that in clude the vector 



nature of polarized scattering (e.g. Stam et al. ( 2006 ): Buenzli Sz Schmid ( 20091 )) instead of the 
scalar approach taken for speed and simplicity in this work, would also be interesting. 

Following the two-stream source function method (jToon et al. I ll989h . we use two-stream quadra- 
ture (sometimes called discrete ordinates) to solve the radiative transfer equation for the upward 
and downward diffuse fluxes. The two-stream quadrature approach is only one of many different 
two-stream methods to approximate the diffuse fluxes. Other approaches, including the hemispheric 
constant, Eddington, modified Eddi ngton, modified quadra t ure, and del ta-fun ction methods, ar e 



treated in detail in sources such as iMeador fc Weaver! <jl98ch . lL.ioJ dl98oll. and 



Toon et al 



(jl989| ). 



In the form of the 'unified description' of two-stream methods in IMeador Weaver! |l98fj), the 
Gaussian quadrature parameters 7« used in this work are: 



(9) 



72 



'3a; 



i-g) 



(10) 
(11) 



74 



1 



73 = 2 (l + VSgnoj 



(12) 



The forward-scattering nature of many atmospheric particulates, especially clouds, is not well 
captured by the two-stream approach as described above. The delta-function technique takes 
advantage of the similarity principle to split out a fractional forward scattering contribution / 
from original scattering phase function into a delta function, and adjust the asymmetry factor g, 
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single scatter i ng al bedo u>, and optical depth r accordingly. This concept is treated in detail in 



Joseph et al.l (|1976l ) who use it al ong with an Eddington two-stream approach. We use a delta- 



function approach as described in iLioul (|1980l ) along with the two-stream source function method 
in this work as well. 



3.2.2. Boundary conditions 

We use the recursive Eq. [5] to calculate the level intensities from the 'ground', N to the top 
of the atmosphere. Note that in discussion of boundary conditions it is common to refer to the 
'ground' or 'surface'. For gas giant exoplanets, these terms are not accurate since there is no 
'ground' in the sense that there is one for a terrestrial planet. Instead, we treat the 'surface' as a 
pressure level deep enough that no photons are reflected back out. We set the following boundary 
condition for the ground level: 

I (r N , »i) = R(»i) ^ (13) 

where F^ total (tjv) = ( t n) + F^ dvrect (tjv) is the total downward flux onto the surface, both diffuse 
and direct. We treat the ground layer as a Lambert surface with reflectance R(fMi) = R- In this 
work, we use a 'black' surface, R = 0. However, in the general case where R ^ 0, we would need 
to solve the diffuse radiative transfer equation to get (tjv). 

We also set boundary conditions at the top of the atmosphere: there is no diffuse flux coming 
from outside the atmosphere, and any upward flux consists of the reflected direct flux off the 
Lambert surface plus the downward diffuse flux scattered upward by the Lambert surface: 

tA — n an d pt — R ( F^ direct + F^ \ (14) 

r top — u clIiu r Surface ~ 1X1 \ r Surf ace ^ r Surface ) \ L ^> 

These are the boundary conditions used with the two-stream source function method described 
above. 



3.2.3. Scattering 



The albedo spectra model uses a wavelength-dependent approximation of anisotrop ic scatter- 
currently 



ing. The cloud model, which uses a full Mie scattering model (jAckerman &; Marley 



2001 



returns the single scattering albedo Co and asymmetry parameter g as a function of wavelength for 
each condensate species and at each temperature and pressure level. In future work, we will update 
the cloud model to return the angular distribution of the scattering in addition to the asymme- 
try parameter. With this type of approach we would then use an M-term fit to the full angular 
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distribution of the scattering function returned by the cloud model in conjunction with a 5— M- 
stream radiative transfer approach. For the moment, we start with the available u and g (which 
are functions of wavelength, temperature, and pressure) as input parameters to a two-term Henyey 
Greenstein function as an approximation to the anisotropic scattering (as discussed below, the 
two-term HG function approximates the angular dependence of both forward and backward scat- 
tering). As with many previous studies, for our current goal of exploring a wide range of exoplanet 
atmosphere models that may be observed by generally photon-starved direct imaging observations, 
our treatment of atmospheric scattering is of a scalar, not vector, nature, meaning that the effects 
of polarization are not considered. We refer readers interested in the effec t s of p o larization to simu- 
lation s of Jupiter-li k e and Rayleigh scattering atmos pheres in IStam et al.l (|2004l ) ; IStam Hovenier 
tod ): IStam et al.l feood ): buenzli fc Schmidl feoogl ). 



We also approximate Raman and Rayle igh scattering . The approximate treatment of Raman 
scattering is the same as that as described in lPollack et al.l ([19861 ). incorporating Ram an scattering 



i nto t he single scattering albedo as a function of the stellar spectrum (see Eq. 2 in lPollack et al 



()1986l )). We approximate Rayleigh scattering by including an additional term in the scattering 
phase function as briefly described below. 

The scattering phase function p (cos 0) is a representative function that takes contributions 
from a number of different scatterers in the atmosphere into account. The scattering phase functions 
used in this model attempt to represent the cumulative scattering effect of all of the constituents 
in each layer of the model. As noted earlier, the single scattering albedo u> and asymmetry factor 
g used in the scattering phase function are updated by the cloud model and radiative-convective 
model at each layer and for each wavelength with weighted contributions from the atmospheric 
constituents. 

The scattering angle can also be written in terms radiation incoming from (p,',(j)') that is 
scattered into (//, <f>): 



cos = fXfjf + (l 



ft 



2\V2 



1 



COS 



(15) 



The second term on the right hand side is sometimes neglected in the case of azimuthal 
independence or averaging. The scattering phase function may be written in terms of cos or 
/j! ', <f>' and /i, <p. The integral of the scattering phase function is typically normalized either to unity 
or to the single scattering albedo, Cj. Note that the scattering phase function p should not be 
confused with planet phase function $ (A, a). Also note the ir difference in directionality between 
the scattering angle and the planet's phase angle, a = tt — 0; backscattering would be more strongly 
seen by an observer at full phase, a = 0°, whereas backscattered = 180°. 

As shown in Fig. [TJ for these models, the contribution of scattering to the albedo spectra not 
only depends on the mixture of gases in a given layer, but also the presence, composition, location, 
and depth of clouds in the atmosphere. There are many different analytic scattering phase functions 
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that attempt to realistically capture the scattering behavior of atmospheric constituents, such as 
the Henyey-Greenstein function, as well as approaches to capturing strong forward scattering using 
delta functions. While there are a variety of scattering phase functions that are used to represent 
different types of atmospheric constituents, it is common to expand the phase function in terms of 
Legendre polynomials in cos0. 



N 



p (cos 8) = "iPt ( cos ) ( 16 ) 



£=0 



21 + 1 f 1 

u t = —I— p (cos 9) Pi (cos 6) (17) 



Note that despite the similar notation, the single scattering albedo, Si, is different from the Legendre 
constants u>e- In most two-stream models, the Legendre polynomial expansion of the scattering 
phase function is not used beyond the first order, Pi. In that case Po (cos0) = 1 and Pi (cos0) = 
cos 0. For a second order expansion, Pi (cos 0) = \ (3 cos 2 — l) . 

For a first order expansion, azimuthal independence is often assumed, leading to: 



p{n,n') =l + 3gnn' (18) 

where g is the asymmetry factor, defined as the first moment of the phase function, and also equal 
to u)i/3. The asymmetry factor g represents the relative strength of forward scattering. 

The Rayleigh scattering phase function is: 

p (cos 0) = - (1 + cos 2 0) (19) 

The cos 2 dependence is not included in the first-order expansion in Eq. [THl However, continuing 
the Legendre expansion to second order allows representation of Rayleigh scattering with the cos 2 
term: 

p (cos 0) = Cjq + Six cos + y (3 cos 2 0- 1) (20) 

For both isotropic scattering and Rayleigh scattering, g = 0. We use the original first order 
expansion (Eq. ITS]) for p in the initial radiative transfer calculation of the two-stream fluxes that 
are to be used in the source function method; the truncation treats the Rayleigh scattering as 
isotropic (since for Rayleigh, g = 0). However, once we have used the source- function method to 
plug the fluxes back in, we use a second-order phase function expansion for p and assume azimuthal 
symmetry such that: 
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p (/if/) = 1 + 3 W ' + ^ [3fi 2 [(if) (21) 

With this expansion, Rayleigh scattering can be modeled when g = and gi = \- In the model, g 
goes to zero when the contribution from TR ay i e i g h in the denominator is large, and #2 goes to 1/2. 
For the purpose of evaluating the second-order expansion with the two-stream source function, we 
also define a 'mean square angle' p,2'- 

2 _ f(l 2 Id(l 
^~ fldfi 

where JI2 is treated as a free parameter ranging from to 1 that depends on the two-stream 
angle approximation used, and that, in combination with the higher order term g2, is intended to 

1 /2 

provide an accurate result for the pure Rayleigh case. In the model used, /Z2 ~ (l/y/oj . The 
weighted asymmetry parameters (5,52) i n the model balance the relative contribution of Rayleigh 
scatterers such that the diffuse scattering phase function approaches a Rayleigh scattering phase 
function when Rayleigh scattering dominates. This capability becomes more important at short 
wavelengths and in clear (cloud and haze- free) atmospheres, and is particularly relevant for the 
direct imaging of exoplanets in the optical. 



.2 r..f\ 2 



3.2.4- Two term Henyey Greenstein 

We use a two term Henyey-Greenstein function (TTHG) as the direct scattering phase function 
to capture the contributions of very strong forward scattering and moderate back scattering. The 
TTHG makes use of the asymmetry parameter g that is calculated in the radiative -convective 
model for the clouds that form within each layer following lAckerman &: Mar ley! (120011 ). As noted 
above, we plan to update this approach to retain the angular properties from each layer and use an 
M-term Legendre polynomial to best fit the angular distribution of the scattering phase function 
for each layer. For the moment, we use a TTHG that is a function of g for each layer and that has a 
form generally consistent with high forward scattering and moderate back scattering by clouds and 



aerosols, a s shown in Fig . 



planets in iPollack et al 



5] and informed by studies of the scattering function in solar system giant 



19861) ; iDvudina et al.l (|2005l ) and discussion of the scattering properties 



of clouds and aerosols in iLioul (|l980l ). The values of g from the cloud model used in this work, at 
pressure levels where clouds have formed, typically range between 0.5 at short wavelengths to 0.99 
at long wavelengths. 

The general form of a single Henyey-Greenstein function is: 



i-g 1 



(1 + g 2 - 2 5 cosG) 



m 3 /2 



(23) 
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The HG phase function is useful because it both approximates strongly peaked forward scattering 
and has a simple expansion in Legendre polynomials: 



N 



Prg = £(2£ + l)<M(cose) 



(24) 



1=0 



From its expansion, pug has a first moment of g, and has a second moment / = g 2 . The ability 
to expand the HG function in terms of N Legendre polynomials as well as the simple form of its 
second moment also makes it useful for adjusting two-strea m solutions to bett er include the strong 
forward scattering nature of cloud particles and aerosols (jJoseph et al.lll976l ). However, a single 
term HG scattering phase function does not capture backscattering. To capture backscattering, 
two term HG scattering phase functions can be used: 



Ptthg = frPHG (9a, 6) + (1 - b) p HG (g b , 6) 



(25) 



where g a > causes the forward peaks and < the backward peaks. The form used here has 
9a = 9, 9b = -g/2, and b = 1 - g\. 



3.2.5. Planet Phase Geometry 

To compute the albedo sp ectra, we integrate over emergent inten sities resulting from the 
incident radiation. Previously ()McKay et al.lll989l : iMarley et al.1 Il999al ). the emergent intensity 
calculations assumed an a = 0° system geometry where the source and observer were collinear, 
such that fj,Q = fix. With that restriction, to consider variation with phase, the geometric albedo 
spectra could then be multiplied by some planet phase function <3> (A, a) to emulate the change in 
emergent flux with phase. Such an approach does not capture the changes in the albedo spectra with 
wavelength that occur for different phases as the paths through the atmosphere and contributions 
from anisotropic scatterers change. We thus have updated the albedo spectra model to use different 
/xo and 

To simulate a spherical planet, we cover the illuminated surface of a sphere with many plane- 
parallel facets, where each f acet ha s differ e nt incident and obser ved angles, as shown in Figure 
[6l Following the approach of iHorakl (119501 ) ; iHorak &; Little) (119651 ) , we use two-dimensional plan- 
etary coordinates and Chebyshev-Gauss integration to integrate over the emergent intensities and 
calculate the albedo spectra. The Chebyshev-Gauss angles are the input angles to the intensity 
calculations. The plane-parallel approac h currently does n ot accurately capture curved geometry 
(limb effects) at high phase angles, e.g. Chapman ( 1931al lbh : we will include this refinement in a 
future update. 
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4. Results 

In this section, we present the albedo spectra model and place our results in context with 
previous work. We first examine the cloud structures as a function of planet-star separation and 
consider their effect on the optical depth r. Next, we present the albedo spectra for all of the 
exoplanet models at full phase (a = 0°). Then we examine how the albedo spectra change as a 
varies from 0° to 180° for representative cases. 



4.1. Clouds and Planet-Star Separation 

In Figures [7J and El we show how the optical depth r varies throughout the model atmospheres 
with increasing planet-star separation for our 3x Jupiter and 10 x Neptune models. The corre- 
sponding albedo spectra are shown in Fig. [9l Figures [7J and ?? shows pressure P vs. wavelength 
with shaded contours of logr that correspond to the color scale on the right. Optical depth of 
r = 1 is at logr = 0. 

For the cloud-free 3x Jupiter case at 0.8 AU, Fig. [7J shows that Rayleigh scattering dominates 
the short wavelengths deep into the atmosphere. Contributions from the alkalis, Na (~ 0.59 //m) 
and K (~ 0.78 //m, see Table [3]) shape the spectrum and are evident deeper in the atmosphere 
where they are favored by chemical equilibrium. Thus even though the effective temperature of 
this model (Table [1]) is 290 K, well below the temperature at which we would expect the atmosphere 
to be devoid of alkali elements, alkalis remain an important opacity source. The reason for this is 
that the clear atmosphere allows photons to penetrate as deep as 10 bars before they are scattered 
for the 0.8 AU 3x case (Figure top left panel). At ten bars the temperature exceeds 800 K and 
the alkalis are present in the gas and thus are detectable. 

At 2 AU, bright and optically thick water clouds form below ~ 500 mbar; their presence 
dominates opacity across the optical. Mie scattering predicts that the extinction efficiency at 
optical wavelengths will be essentially constant for all particles with radii greater than 1 /im and 
thus the spectral contribution of the clouds appears flat in the figures. It is interesting to note 
that optical thicknesses can reach ~ 10, 000 for H2O clouds. At 5 AU, ammonia clouds form below 
~ 50 mbar and water clouds form below ~ 1 bar, at a deeper level than the for the warmer 2 
AU case. The ammonia clouds also contribute significantly to the optical depth. For the still 
colder 10 AU Jupiter 3x model, both the NH3 and H2O clouds form deeper into the atmosphere, 
allowing Rayleigh scattering to again dominate the opacities at short wavelengths. Even at long 
wavelengths, unit optical depth is likely to be reached before encountering the clouds. 

In addition to the alkalis (clearly seen at 0.8 AU) and CH4, opacity due to gaseous H2O plays 
a role between ~ 0.92 and 0.95 /im in our albedo spectra for all of the cases closer in than 5 AU. 
Detection of water absorption features in the optical have not been con firmed in Jupiter's reflection 



spectrum, although absorption features near ~ 0.94 fim were noted by iKarkoschkal ()1994! ) as being 
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present. At the time, H2O was suggeste d potentially bein g the cause, however, later observations 



suggested that NH3 could be responsible (|Karkoschkalll998l ). We address the presence of H2O again 
later in discussion of Figure [TT1 

Figure Oshows the relationship between r, P, A and planet-star separation for the 10 x Neptune 
models. Similar to the 0.8 AU 3x Jupiter case, the 0.8 AU 10 x Neptune case is dominated by 
Rayleigh scattering at short wavelengths, only with more pronounced absorption features into the 
red. The 2 AU case is also dominated by the presence of H2O clouds that form just below ~ 380 
mbar, a bit higher than the Jupiter case, implying that the albedo spectra should be a bit brighter 
as well. At 5 AU, an ammonia cloud forms below ~ 100 mbar, and extends down nearly to the 
H2O cloud, which forms below ~ 470 mbar. Also similar to the 10 AU 3x Jupiter case, the 10 AU 
10 x Neptune shows both cloud decks forming at substantially higher pressures, below ~ 700 mbar 
for NH3 and below ~ 3.6 bar for H2O. 



4.2. Geometric Albedo Spectra (a = 0°) 

In Figure [9] we show albedo spectra for all of the model cases summarized in Table [H The 
spectra f or our Jupiter-li k e exop lanet models can be c o mpar ed with those presented in lMarlev et al 



rapi 

(ll999aT ): ISudarskv et al.1 (120001 ) and ISudarskv et all (12005T ). For the aurpose of comparison, the 



model Jupiters used here map to the Class III (clear), Class II (water cloud) and Class I (ammonia 
cloud) nomenclat ure used in ISudarskv et al.l (|2000l ). Our clear and ammonia cloud models are 
similar to those in ISudarskv et al.l (J2000), and our water cloud models at 2 AU are a bit brighter. 
As shown in Fig. [lOj distinct Na, K, CH4 and H2O features are apparent, particularly CH4 near 
0.62, 0.74, and 0.89 /mi (see Table [3]). As discussed in p~T| the cloud-free 0.8 AU spectra are 
dominated by Rayleigh scattering at short wavelengths for both lx and 3x Jupiters and 10 x and 
30 x Neptunes. 

At any given planet-star separation, the higher metallicity compositions generally have smaller 
albedos due to the increased opacity of their atmospheres. The presence of relatively high and thick 
H2O clouds in the atmosphere of the 2 AU Jupiters result in a higher albedo across the visible. The 
bright effect of relatively high H2O clouds is also apparent through mid-band in the 2 AU Neptunes, 
however, the absorption features are considerably more pronounced at longer wavelengths for the 
Neptunes than for the Jupiters. At full phase, backscattering from clouds also plays a role. At 
separations of 5 AU and 10 AU, the albedos show progressively larger contributions by Rayleigh 
scattering at short wavelengths, consistent with the lower clouds in Figs. [7]and[8l At 5 AU and 10 
AU, the difference between albedos for lx and 3x solar abundances of heavy elements for Jupiters 
is larger at all wavelengths than the difference between albedos for 10 x and 30 x solar for Neptunes. 
In the cooler atmospheres at larger planet-star separations the albedo spectra seem to change little 
with increasing heavy element abundances above about a ten-fold enhancement. 



We compare the albedo spectra for our standard 5 AU 3x enhancement Jupiter model at a = 0° 
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with o bserved data from the real Jupiter in our solar system at near full-phase from iKarkoschka 
(| 19941 ) in Fig. [TTJ This figure illustrates how our interpretation of Jupiter might proceed if we 
were to detect it as an exoplanet. The general agreement in morphology of the spectral features 
is simply a consequence of the spectrum being primarily (but not exclusively) shaped by methane. 
As noted in Q3.2\ we do not include the effect of photochemical products such as haz es that would 



explain the difference betwee n our model and the observed data at short wavelengths (|Marley et al 



1999a 



Sudarsky et alJ I200Q ). Although our 5 AU 3x Jupiter model was not adjusted to fit the 



data, varying the cloud thickness or / se( j parameter would brighten or darken the spectra. 



The difference between our 5 AU 3x Jup iter model and th e dat a from iKa rkoschkal (|1994l ) 
around 0.94 /im are of interest. As noted above IKarkoschkal (|1994l ) and Karkoschka (1998) noticed 
features near this wavelength region in Jupiter's albedo spectrum. IKarkoschkal (| 19941 ) mentioned 
that they could potentially be water features, but in the later paper concluded they were more 
likely to be ammonia. W hile the features w e see in this regi o n of o ur models are likely not the 
same as those observed by IKarkoschkal (| 19941 ) and IKarkoschkal (119981 ) . we were able to confirm, by 
running models both with and without H2O as shown in Figlll[ that the model features around 0.94 
/^m are indeed attributable to H2O. These features would be interesting targets for direct imaging 
observations of exoplanets since the appearance of water is a sensitive probe of temperature and 
composition. 



4.3. Albedo Spectra vs. Phase 

In Fig. [12] we show the albedo spectra of the 0.8 AU 3x Jupiter model as a progresses from 0° 
(full) to 180° (new). As expected, the albedo decreases with increasing a. More importantly, Fig. [12] 
also shows the ratio of the albedo spectra at each a increment to the geometric albedo spectrum 
at a = 0°. The wavelength dependent changes are apparent as a function of a. There are both 
broad effects on the spectrum, such as the drop in brightness from shorter to longer wavelengths 
becoming shallower as the planet phase increases, as well as narrower effects on the spectrum due 
to the changing line depths as different volumes of atmosphere are probed at different phases. 

Figure [TBI shows albedo spectra of the Jupiter and Neptune models at 5 AU for a few different 
phase angles to illustrate the point noted in §4.21 that for planet-star separations beyond 2 AU, 
where it is cooler, there is not as a large of a difference in the albedo spectra for the 10 x and 30 x 
Neptune models as there is between the lx and 3x Jupiter models, even at different phase angles. 
This implies that it will be more challenging to use albedo spectra to differentiate between Neptune 
analogs with large abundances of heavy elements. In future work, we plan to further increase the 
metallicity of Jupiter analogs to determine to what extent this occurs for the Jupiters as well as 
Neptunes. We also plan to investigate decreasing the metallicity of the Neptune analogs to examine 
the range of metallicities that we can differentiate between at cooler, larger planet-star separations. 
For the Jupiter case, it appears that the difference in metallicity is easier to detect when probing 
deeper into the atmosphere at smaller phase angles than at larger phase angles. 
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4.4. Phase Functions 

From the albedo spectra as a function of phase, we can also generate phase functions, <E> (A, a) 
as in Eq. [3j Although we can calculate a phase function for each wavelength individually, we 
instead present the phase functions for our mo del exoplanets using the Johnson-Morgan/Cousins 
UBVRI filter passbands (IFugukita et al.l 11995 ) as shown at the top of Fig. [TH Also shown for 
reference in Fig. [14] with the filter responses^ is the albedo spectrum of the 0.8 AU lx Jupiter 
model. Note that our models extend shortward only to 0.35 ^m. Since they do not extend to the 
shortest wavelengths of the U filter, we do not use U in this work. Fig. [TH also shows the solar 
spectrum n Se d fa the parent stafi and the reference spectrunfl used to compute color magnitudes. 
Direct imaging observations will obviously be made with different filter responses, possibly tuned 
specifically to achieve the particular scientific objectives of the observation. However the approach 
presented here would be applicable to any arbitrary filter set as well (see Appendix IA.1|) . 

In Figure \TE\ we first compare the phase function for our 10 AU 1 x Neptune mode l with 
data points of both Uranus and Neptune from Voyager 1 presented in IPollack et al.l (1986). The 
data points and our phase function are shown with a Lambert phase function for reference. While 
the error bars are relatively large, and there are two Uranus data points that are closer to the 
Lambert curve than our model, there is general agreement at higher phases. Since these particular 
Voyager 1 data were taken with a clear filteilfl, the model phase function shown is an average of 
phase function s over all wavelengths. We did n ot include the phase function for Jupiter used in 
Dvudina et al.l (|2005l ) and lSudarskv et al.l (|2005l ) for comparison since the Pioneer data used in the 
generation of that curve were not disk-integrated. 

Figure [TBI shows phase functions for each model exoplanet and for the B, V, R, and I filters. 
The phase function is shaped by geometry and the composite of all of the scattering phase functions 
making contributions to scattering in the exoplanet's atmosphere. For reference, the area beneath a 
Lambert phase function, or its phase integral q, as in Eq. [H is 3/2, and the area beneath a Rayleigh 
phase function is 4/3 ((/Lambert > ^Rayleigh)- Different line styles represent different metallicities: 
lx is solid, 3x is dashed, 10 x is dot-dash, and 30 x is dotted. A Lambert phase function is shown 
in each sub-plot for reference. For the 0.8 AU cloud-free case, it is difficult to distinguish between 
the metallicities, although the higher metallicity phase functions are generally larger than the lower 
metallicities. The phase functions start off quite similar to the Lambert phase function at short 
wavelengths, and become shallower and rise above than the Lambert phase function as the colors 
move red ward. 



3 Filter responses obtained from the Virtual Observatory, http://voservices.net/filter/filterfindadv.aspx 
September 2009 

4 Solar spectrum obtained from STScI, ftp://ftp.stsci.edu/cdbs/calspec, September 2009 
5 Reference spectrum obtained from STScI, ftp://ftp.stsci.edu/cdbs/calspec, September 2009 
6 The Voyager vidicon detector was sensitive from 0.28 to 0.64 /im. 
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For the 2 AU water cloud case, the phase functions are all just below the Lambert phase 
function; the difference becomes more pronounced as colors move redward; for the R and I filters, 
the beginnings of a forward-scattering 'toe' are seen at large phase angles. Effects of backscattering 
would appear at small phase angles as well. For the 5 AU ammonia and water cloud case, the 
ammonia clouds are higher in the atmosphere than the water clouds, and there is a noticeable 
difference between the Jupiters and Neptunes, likely because the ammonia clouds are a bit higher 
for the Jupiters than they are for the Neptunes (also true at 10 AU, see Figs. [7]and[8]). At shorter 
wavelengths at 5 AU, the Jupiters fall below the Lambert curve and the Neptunes land closer 
to it, although the Neptune curves also fall below the Lambert curve as colors get redder. The 
forward-scattering 'toe' becomes more pronounced. At 10 AU, both the ammonia and water clouds 
condense at much higher pressures (lower altitudes), particularly for the Neptunes. The Neptune 
phase functions are below both the Jupiter phase functions and the Lambert pha se function. We 



note t hat our phase functions are not nearly as steep near full phase as are those in lSudarskv et al 



(|2005l ). We will further examine the dependence of phase function morphology on the assumptions 
made in the treatment of the direct and diffuse scattering functions in a future paper. We discuss in 
the following section the use of color-color diagrams to help distinguish between the different types 
of exoplanet models, as it would be difficult to observe and confirm the small differences between 
the phase functions shown in 1161 

From the phase functions and the geometric albedos, we can use the equations outlined in 
§2.21 to calculate the phase integral q, spherical and Bond albedos A s and Ab- Our Bond albedo 
calculation integrates over 0.35 to 2.5/xm. Contributions from longer wavelengths are negligible 
(less than 5%) for all cases since both the albedo spectra and incident flux drop off quickly at 
longer wavelengths. Table [4] contains geometric albedos calculated for each B, V, R, I color filter 
for each of our models. Table [5] contains the corresponding phase integrals q. The Bond albedos in 
Table [6] yield the effective temperatures in Table CD 



4-4- 1- Bond albedos and implications for dynamics 

The Bond albedos shown in Table [6] combined with our assumed internal heat flows result 
in planetary effec t ive te mperatures for our model planets shown in Table [T] (using equation [7]). 



Ingersoll Sz Pored (|1978l ) found that if the ratio E of emitted to total absorbed power is greater 
than 4/7r ~ 1.27 for a giant planet, then internal convection will redistribute heat such that there 
is vanishing equator to pole temperature gradient, as there is on Jupiter. Based on our computed 
Bond albedos and internal energy fluxes (parameterzed as Tj n t in Table 1) then all of our model 
cases-except for the Jupiters at 5 and 10 AU-are in fact characterized by E < 1.27. Thus we expect 
that in the general case of mature Jupiters and Neptunes (with modest Ti nt ) many planets will fall 
into a regime in which we expect to find substantial equator to pole temperature gradients. Such 
planets may be expected to exhibit atmospheric dynamics and cloud patterns quite unlike those 
expressed by Jupiter and Saturn. 
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5. Application to Observations 



First generation space-based coronagraphic direct imaging observations of widely-separated 
(cooler) exoplanets will likely consist of photometric measurements using fairly wide filters, and 
result in coarse spectra or possibly even only few photometric data points in the visible. Here, 
we use the high-resolution model albedo spectra to generate coarse (R = 5 and 15) spectra and 
color-color diagrams of our gas and i ce gia nt planet models. We compare our model coarse spectra 
and exoplanet colors to iKarkoschkal ()1994l )'s observations of the giant planets in our solar system. 



5.1. Coarse Spectra 

We noted in ^2] that ey e n for the low-resolution spectra of solar system outer planets in Fig. [2] 



generated from IKarkoschkal (|1994l ) 's high resolution spectra, it was still possible to discern between 
gas giant (Jupiter, Saturn) and ice giant (Uranus, Neptune) planets. This is still the case in Fig. [T71 
where the Jupiter models are generally brighter than the Neptunes. For solar system outer planets, 
the real Jupiter and Saturn are darker at short wavelengths due to absorption by photochemical 
hazes, but are brighter in the middle of the optical. The real Neptune and Uranus are brighter 
at short wavelengths and darker at longer wavelengths, due to stronger absorption features in the 
red. The brightness of a given planet depends sensitively on the height and optical thickness of 
the uppermost cloud decks, however, so one could imagine planets with somewhat different cloud 
structures from those presented here. The actual flux measured by a camera depends on the 
product of a planet's albedo, cross section, and phase. Within the limitations of the models (which 
do not include photochemistry) ice giants generally will be more difficult to detect than gas giants 
at longer wavelengths both because of their lower albedos and their smaller radii. 

Fig. [17] again illustrates the difference between lx and 3x solar abundances for the Jupiter 
models is larger than the difference between 10 x and 30 x for the Neptune models. In the R = 15 
spectra in Fig. \T7\ it is also possible to see the effect of the strong methane absorption in all of 
the spectra between 0.7 and 0.9 //m. Placement of relatively coarse filters both centered on these 
features as well as on the continuum around them should make it possible to detect and possibly 
differentiate between these features in early direct imaging observations. 



5.2. Color-Color 



As mentioned in §4,44 to generate illustrative c olor-color diagrams o f exoplanets, we use the 
Johnson-Morgan/Cousins UBVRI filter passbands (IFugukita et al.lll995l ) as shown a t the top of 
Fig. 1141 In Fig. \18\ we first use the high-resolution solar system outer planet data from IKarkoschkal 
Jl99j" ;o compute the colors for Jupiter, Neptune, Saturn, Uranus, and the moon Titan (see 
§A.1I for methods). These serve as a useful reference for comparison with the model Jupiters and 
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Neptunes. We then compare these with the (V — R) and (R — I) colors of the exoplanet models. 
Different symbols represent different planet-star separations: circle for 0.8 AU, square for 2 AU, 
triangle for 5 AU and diamond for 10 AU. Red represents the lx Jupiters and magenta the 3x 
Jupiters. Blue represents the 10 x Neptunes and cyan the 30 x Neptunes. The symbols are strongest 
at a = 0° and fade with a increasing to 180° in 10° increments. 

The cloud- f ree l x model at 0.8 AU exhibits behavior that is generally consistent with that in 



Sudarskv et all (|2005l ). ranging from (R - 1) of about -1 to -0.5 and (V - R) of about -0.4 to 0.2 
with increasing phase. The 0.8 AU 3x Jupiter model and both Neptune models are also consistent 
with this behavior. The lx and 3x Jupiter models at 2 AU and 5 AU at full phase are quite close 
to the real Jupiter. The 10 x and 30 x Neptune models also approach the observed Neptune, which 
is encouraging given that the brightness for the models are a bit darker than the observational 
data. We also note that for the cloudy Jupiter models, it is difficult to differentiate between the 
models at planet-star separations wider than 2 AU on the (V — R), (R — I) color-color diagrams. 

While additional work remains to be done on both generating and interpreting exoplanet albedo 
spectra models, another relevant question for color-color comparisons is which colors (parts of the 
spectrum) are most informative. Making color-color comparisons using the filters that best allow 
differentiation between planet models also requires further optimization. In Fig. [19] we show four 
additional color-color diagrams using different filter pairings. From this figure, it appears that the 
(B — V), (R — I) and (B — V), (V — I) diagrams provide the greatest spread between the planet 
types. These diagrams are also useful for determining how well the models capture the integrated 
effect of different broad regions of the spectrum in comparison to the solar system outer planet 
data. 



It has been suggested by iTraubl (|2003l ) based on the upper left panel of Figure [18] that color 
alone will be a definitive test of planet characteristics, for example by distinguishing terrestrial, 
ice-giant, and gas- giant planets. Figures [T8l and [T9l suggest that such an approach is not as straight- 
forward as it might appear. Because planetary phase functions vary with wavelength the planet 
colors also depend upon phase. Furthermore, colors vary as cloud heights vary with atmospheric 
temperature and composition. Thus for the example at 5 AU in the (V — R) vs. (R — I) color space, 
the distinction between the Neptune and Jupiter-like models is slight, particularly if the phase is 
unknown. We conclude that a single color-color diagram is likely inadequate to sort planets with 
certainty and the more photometric bandpasses that are available, the better. Of course, the best 
characterization of an exoplanet will be with spectral information, preferably with i? > 15 (Figure 

QZD- 
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6. Conclusions and Future Work 
6.1. Conclusions 

In this work, we have investigated the roles of planet-star separation, heavy element abundance, 
and planet phase on the reflected-light albedo spectra of Jupiter and Neptune analog exoplanets. 
Albedo spectra of gas giant exoplanets are sensitive to both the presence of clouds as well as to 
the changes in optical depth that occur as a function of planet phase, when different volumes of 
atmosphere are sampled corresponding to the source-observer geometry. The planet-star separation 
(temperature) plays a primary role in determining the cloud structure in our model exoplanet 
atmospheres. The abundance of heavy elements (metallicity) plays a supporting role, with clouds 
typically extending deeper to higher pressures and higher opacities with increasing metallicity. 



Mar lev et al.l (|1999al ) , ISudarskv et al.l (|2000l ) , and ISudarskv et al.l (|2005l ) have previously ex 



plored the effect of changing separation on albedo spectra of 1 Mj and more massive planets. We 
expand upon this work, developing a different model implementation to include observational geom- 
etry, exploring the effects of changing atmospheric composition and adding Neptune-mass planets. 
In total, we present sixteen different model cases that span a range of parameter space consistent 
with the expected capabilities of early direct imaging observations of exoplanets in reflected light. 
The sixteen cases include a range of planet sizes (Jupiter and Neptune analogs), planet-star sepa- 
rations (0.8 to 10 AU), compositions (heavy element abundances between lx and 30 x solar), and 
planet phases (from 0° to 180° in 10° increments). The exoplanet models capture cloud-free, H2O, 
and H2O + NH3 cloud cases. 

As previously established, clouds play a key role in controlling the brightness of a planet seen 
in scattered light. Water clouds first form high in the atmosphere for planets found beyond about 
1AU. With increasing planet-star separation, the water clouds are found progressively deeper in 
the planetary atmosphere until ammonia clouds form above them at planet-star separations ~ 5 
AU. For cloudless planets, gaseous alkali elements, particularly Na and K, are important opacity 
sources along with water. Once clouds form, the main absorbing species is methane. However, the 
reflected light spectrum of even Jupiter itself is not purely a me thane spectrum, for example, it also 



contains ammonia absorption features (|Karkoschkalll994l . Il998l ). Our computed energy budgets for 



our model mature gas and ice giant planets also imply that there should be substantial equator 
to pole temperature gradients and thus atmospheric dynamics that will differ from solar system 
giants. 

We find that, as in the solar system, the reflected planetary spectra are sensitive to planetary 
phase as planets are not simply Lambertian scatterers at all wavelengths. As a result, knowledge of 
the phase at which the image of an exoplanet was obtained is a key ingredient to the interpretation 
of the observation. Planet phase will be trivially known for planets that have also been detected 
by radial velocity observations. It will be more difficult to interpret observations of exoplanets seen 
in only one or two images (perhaps because the planet has passed inside or outside of the inner or 
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outer working angles of the coronagraph) and for which no RV data is available. 



Ratios of albedo spectra at taken different planet phases clearly have wavelength-dependent 
features. Making multiple measurements of exoplanets at different phase angles, taking the ratio, 
and examining the relative changes in the depth of absorption features can help to characterize the 
vertical distribution or scattering dependences of atmospheric constituents in an exoplanet's atmo- 
sphere. The spectra generated for these exoplanet models can be used to estimate the sensitivity 
required to detect this type of phase-dependent effect in direct imaging observations. 



As has been noted elsewhere (e.g.. lMarleyill997l : iMarley et al.ll2007l ). photochemical hazes will 
add a challenging aspect to any effort to interpret cool giant exoplanet spectra. As the comparison 
of Jupiter and the model Jupiter in Figure HU attests, the neglect of photochemical hazes is a first 
order discrepancy in the model. A Neptune at 1 AU around a solar type star, with a rich brew of 
not only methane, but also water and H 2S in its atmosphe re will unquestionably exhibit a complex 
photochemistry and thermochemistry. IZahnle et al.l (|2010l ) has only begun to explore the chemical 
processes that might take place in such atmospheres, but regardless we can expect that actual 
spectra may deviate strongly from those shown here. If so this will add new layers of insight to 
that gained from the study of the scattering spectra of giant exoplanets. 

Coarse spectra and color-color diagrams are a first step in characterizing an exoplanet's atmo- 
sphere. We use standard UBVRI filters to calculate the colors for each of our exoplanet models at 
different phase angles as well as for solar system planets Jupiter, Saturn, Saturn's satellite Titan, 
Neptune, and Uranus. Exoplanet colors are highly dependent on planet-star separation (clouds) 
and planet phase, and less sensitive to metallicity, although it may be possible to differentiate 
between high and low metallicity for certain filter pairings. It is important to make multiple obser- 
vations and to obtain as much information as possible about the geometry of the system in order 
to interpret the observed spectra. 



6.2. Future Work 

There are several areas that could be further developed based on this initial work. The models 
and methods presented here can be used to study systems using different stellar spectra and with 
different instrument-specific filters. Signal to noise ratio estimates can then be computed for spectral 
features of interest, such as Na, K, CH4, NH3, and H2O at appropriate resolutions to help inform 
instrument specifications and observing plans. The exoplanet albedo spectra can be combined with 
model spectra of exozodiacal dust and used as input to instrument performance simulations to 
investigate the detectability of exoplanets as a function of wavelength in systems with substantial 
circumstellar debris structure. 

Additional exoplanet model 'points' in the metallicity and planet-star separation space can 
be considered. Given the relatively small effect of increasing from 10 x to 30 x solar metallicity 
on the albedo spectra for Neptunes at separations larger than 2 AU, and yet the noticeable effect 
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of increasing from lx to 3x solar metallicity for the Jupiters, it would be interesting to adjust 
abundances and planet-star separations to determine where the effect begins to diminish. 

Cloud formation has a strong effect on albedo spectra. It would be instructive to consider the 
effect of different sedimentation ratios for these exoplanet models. It would also be interesting to 
adjust the vertical extent of the clouds to investigate the effect of longer path lengths through the 
clouds at different phase angles and for different vertical distributions. A similar investigation to 
compare the absorption features of well-mixed versus stratified or pressure-dependent atmospheric 
species as a function of phase would also be useful in preparing to characterize exoplanets with 
direct imaging observations. These types of investigations would benefit from updating the current 
cloud model to output the full anisotropic phase function in addition to the asymmetry parameter 
and single scattering albedo. 

In this work, we restricted the system geometry to the simple case of circular orbits with no 
other variation in the Keplerian elements (such as eccentricity, inclination, or obliquity). Investiga- 
tion of more complex orbits and the effect of the changing planet-star separation through the orbit 
on the T(P) structure of the atmosphere would be informative. The time-dependence and evolution 
of an exoplanet's atmosphere through a more complex orbit could have a substantial effect on the 
observed albedo spectra. The obliquity of a planet and the potential for seasonal change would 
further complicate the interpretation of albedo spectra. 

With the ability to create many facets over an illuminated sphere, we could use models for sev- 
eral types of facets to paint a portrait of a non-uniform exoplanet on the sphere. The contributions 
from the different facets would then be integrated to generate the albedo spectra. For example, we 
could create facets that correspond to banded structure, or other variations in cloud topology. 



Relative color-color comparisons, such as u — v , can be made using the filter response, the 
albedo spectrum, and parent star spectrum: 



where Ru,a and Ry,A are the U and V filter responses as shown in Fig. Q31 F &t x is the spectrum of 
the parent star, and A a> ^ is the albedo spectrum of the exoplanet at some phase a. Usually it is 
assumed that a = 0° and the geometric albedo spectrum is used, but in this work we are interested 
in changes with a and use the albedo spectra corresponding to different phases. 

Planets in our solar system have apparent magnitudes that assume the planets are illuminated 
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by the sun and observed from Earth. For the purpose of comparison, we can take the approach 
of calculating magnitudes for the exoplanet models as though they were also in the solar system, 
using their geometric albedo spectrum, the solar spectrum, the filter response, and the reference 
spectrum. For example, to calculate V for an exoplanet, 



where -F re f,A is the reference spectrum, A is the distance from the observer on Earth to the planet, 
r is the distance from the parent star to the planet (both in AU) , R p is the radius of the planet (in 
AU) and r is 1 AU. 

To go from planetary to absolute magnitudes (at 10 pc) requires adding 51og 10 (10 • 206, 265), 
where 1 pc ~ 206, 265 AU. 
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Fig. 1. — Simulated direct images of Jupiters and Neptunes around a solar analog at 10 p c , as a 
function of planet-star separation (1 AU, 2 AU, and 5 AU) and A (350 nm, 550 nm, and 750 nm, 
each with 100 nm bandwidth). The simulation is for a D = 4 m space telescope with a Phase- 
Induced Amplitude Apodization (PIAA) coronagraph (jGuvon et al.ll2005l ). The PIAA coronagraph 
has an inner working angle ~ 2X/D. The integration time is 10 hours and effects from photon noise 
and 1 zodi of both local and exozodiacal dust are included (system inclination is 60°). Note the 
scale with A in the se 256 x 256 pixel i mages: 350 nm ~2.8 mas/pix, 550 nm ~4.4 mas/pix, and 750 
nm ~6.0 mas/pix (jCahoy et al.ll2009l ). 
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Fig. 2. — Low-resol ution albedo spectr a of the outer planets in our solar system. Original high- 
resolution data from iKarkoschkal ( 1994 ). Left, R = A/AA = 5, right, R = 15. 
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Fig. 3. — Temperature and pressure models used as input to the albedo spectra model. Left, models 
for Jupiter with compositions of lx and 3x solar and at separations of 0.8 AU, 2 AU, 5 AU, and 
10 AU. Right, models for Neptune with compositions of 10 x and 30 x solar and at separations of 
0.8 AU, 2 AU, 5 AU, and 10 AU. See additional model descriptions in $3]and Table [TJ 
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Boundary 



Fig. 4. — Diagram of the plane-parallel geometry for the albedo spectra model, where [Iq represents 
the cosine of the angle of incidence, f3, and \x\ represents the cosine of the angle of reflection, 
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Fig. 5. — The two-term Henyey-Greenstein phase function used in this work, shown with respect 
to both the scattering angle as well as phase angle a for reference. Backscattering is at = 180° 
(a = 0°) and forward scattering at = 0° (a = 180°). 
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Fig. 6. — Each dot represents one plane-parallel albedo spectra model run with a different pair of 
incident /io and emergent fi\ angles. A two-dimensional Chebyshev-Gauss integration over all dots 
is performed to calculate the albedo spectra as a function of planet phase. There are n = 100 dots 
shown per planet phase. 
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Fig. 7. — These plots show how the optical depth, r, varies as a function of pressure and wavelength 
for the Jupiter 3x models. The background shading represents logr, with values corresponding to 
the color scale on the right. Optical depth of 1 is at logr = 0. The H2O and NH3 cloud structures 
represent the location of the clouds in the models. See £ )4.1l 
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Fig. 8. — As in Figure [7| but for the Neptune lOx models. 
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Fig. 9. — Geometric albedo spectra (a = 0°) for each of the exoplanet model atmospheres used in 
this work. The models cover a range of planet-star separations from 0.8 AU to 10 AU, and a range 
of heavy-element abundances (metallicities) with respect to solar (lx). The Jupiter models have 
lx (solid) and 3x (dashed) solar abundance, and the Neptune models have 10 x (solid) and 30 x 
(dashed) solar abundance. See §4. II and 14.21 
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Fig. 10. — Geometric albedo spectra of Jupiter-analogs at 0.8 AU (red) and 2 AU (green) and lx 
(solid) and 3x (dashed) solar heavy element abundances; prominent spectral features are noted: 
CH 4 , K, Na, and H 2 0. See Tabled 
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Fig. 11. — Left, comparison of albedo spectra from our 'untuned' 5 AU Jup iter-mass exoplane t 
model (3x solar abundance of heavy elements) with data of Jupiter obtained by iKarkoschkal ()1994l ). 
See §4.21 Right, identification of the features near 0.94 /im in our model: they are water features. 
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Fig. 12. — Albedo spectra of the 0.8 AU 3x Jupiter as a function of planet phase a. Left, albedo 
spectra vs. a. Right, ratio of albedo spectra in increments of a = 10° to that of the albedo 
spectrum at a = 0°. Note the variation in the ratio as a function of wavelength; see £14,31 
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Fig. 13.— The effect of metallicity as a function of a at 5 AU for a = 0°, 30°, 60°, 90°, 120°, and 
150° (blue, red, green, cyan, magenta, and yellow, respectively). Left, Jupiters, solid lines are lx 
solar, and dashed are 3x. Right, Neptunes, solid lines are lOx solar, and dashed are 30x. Note 
that at 5 AU, the difference in heavy element abundance has a larger effect on the albedo spectra 
for Jupiters than it does for Neptunes. The lx Jupiter model is substantially brighter than the 
3x, likely due to the H2O clouds forming at lower pressure (higher altitude). 
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Fig. 14. — Reference spectra and color filters used in the models and for calculating magnitudes and 
for color-color comparisons. Upper, the Johnson- Morgan/Cousins UBVRI filter passbands courtesy 
the Virtual Observatory, shown with albedo spectrum of model Jupiter at 0.8 AU at full phase (black 
line). Middle, the solar spectrum from the Space Telescope Imaging Spectrograph (STIS) courtesy 
the Space Telescope Science Institute (STScI). Bottom, the Vega reference spectrum (also courtesy 
STIS/STScI). See MandlATl 
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Fig. 15. — Comparison of our Neptune 10 AU at 10 x solar metall icity model with da ta at higher 
phase angles from Voyager 1 observations of Uranus and Neptune ([Pollack et al.lll986l ). 
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Fig. 16. — Phase functions shown versus planet-star separation and for different color filters (see 
Fig. [TUfor filters). Each plot shows the Lambert phase function (gray line) as a reference along 
with phase functions for the lx (solid line), 3x (dashed line), 10 x (dash-dot), and 30 x (dotted) 
exoplanet models. 
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Fig. 17. — Simulated coarse-resolution albedo spectra for Jupiters and Neptunes at 5 AU with 
different compositions at a = 0° for comparison with Figure [2j Left, resolution R = 5. Right, 
resolution R = 15. 
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Fig. 18. — Color-color comparisons. Upper left, solar system outer planets colors, calculated from 
Karkoschka (1994). Red = lx and magenta = 3x Jupiters. Blue = 10 x and cyan = 30 x Neptunes. 
Upper right, 0.8 AU vs. a, fading from a = 0° to a = 180° in 10° increments. Middle left 2 AU, 
middle right 5 AU, lower left 10 AU, and lower right, all. 




Fig. 19. — Color-color comparisons, with 'all' shown as in Fig. 1181 lower right, but for comparisons 
between different color filters. The B-V, R-I comparison shows substantial horizontal and vertical 
offsets among the different separations and compositions. 
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Table 1: Description of the model gas and ice giants used as inputs to the albedo spectra simulation. 
The separations and cloud types are noted. T[ nt is a parameterizaton of the flux from the interior 
of the planet, due to its cooling with time (temperature in the absence of any incident flux). f se( i 
is the ratio of the microphysical sedimentation flux to the eddy sedimentation flux. 
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Table 2: Characteri s tic 'e ffective radii' (/zm) of cloud condensate particles predicted by the 



Ackerman Marlevl (|200ll ) model. The parameter 'effective radius' is used for comparative pur- 



poses to summarize the level-dependent distribution of radii actually used in the model 
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Table 3: Approximate wavelengths of optical absorption features noted in our gas giant exoplanet 
models. The references in this table refer to work by other researchers who have noted these 
features. The use of 'weak' refers only to the appearance of the spectral features in this work; for 
example, resolution here is not sufficient to resolve expected doublet absorption features. 



Approximate A (//m) Species Reference 

O40 K 3/7 

0.46 CH 4 1 

0.48 CH 4 1 

0.54 CH 4 1 

0.59 Na 'doublet' 3, 4, 5, 7 

0.62 CH 4 1 

0.65 H 2 weak 3, 6 

0.73 H 2 weak 3, 6 
0.77 K 'doublet' weak 3, 4, 5, 7 

0.78 CH 4 1 

0.79 CH 4 1 

0.83 H 2 weak 2, 3, 6 

0.84 CH 4 1 

0.86 CH 4 1 

0.89 CH 4 1 

0.91 CH 4 1 

0.94 H 2 1, 2, 3, 4, 6 

0.99 CH 4 1 



^arkoschka ( 


1994) 


' Mar lev et al. ( 


L999ah 


' Sudarskv et al. 


(200Ch 


4 Burrows et al. 


(2004) 


* Fortnev et al. ( 


2008b) 



R. Freedman, pers. comm. 

7 NIST atomic spectra database, http://www.nist.gov/physlab/data/asd.cfm 
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Table 4: Geometric albedos (A g ) of the exoplanet model Jupiter and Neptunes in this work, for 
each of the B, V, R, and I filters shown in Fig. [14] and the albedo spectra model described in §3.2.51 
and Appendix [Al 

B V R I 

Model A g A g A g A g 

0.8 AU, Jupiter lx 0.556 0.322 0.156 0.047 
0.8 AU, Jupiter 3x 0.482 0.241 0.102 0.029 
0.8 AU, Neptune lOx 0.455 0.209 0.074 0.016 
0.8 AU, Neptune 30x 0.359 0.142 0.045 0.010 



2 AU, Jupiter lx 0.733 0.742 0.727 0.673 

2 AU, Jupiter 3x 0.758 0.766 0.726 0.630 

2 AU, Neptune 10 x 0.744 0.728 0.616 0.420 

2 AU, Neptune 30x 0.735 0.674 0.498 0.267 

5 AU, Jupiter lx 0.609 0.567 0.531 0.453 

5 AU, Jupiter 3x 0.567 0.506 0.458 0.388 

5 AU, Neptune lOx 0.434 0.326 0.238 0.145 

5 AU, Neptune 30 x 0.430 0.303 0.191 0.090 

10 AU, Jupiter lx 0.450 0.386 0.358 0.313 

10 AU, Jupiter 3x 0.316 0.260 0.246 0.240 

10 AU, Neptune lOx 0.388 0.295 0.228 0.147 

10 AU, Neptune 30x 0.388 0.279 0.189 0.096 
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Table 5: Phase integrals (q) of the exoplanet model Jupiter and Neptunes in this work, using the 

filters shown in Fig. [T4l and the albedo spectra model described in §3.2.51 and Appendix 1X1 

Model g B qv qr qi 

0.8 AU, Jupiter lx 1.470 1.582 1.693 1.764 

0.8 AU, Jupiter 3x 1.501 1.636 1.739 1.784 

0.8 AU, Neptune lOx 1.509 1.639 1.725 1.748 

0.8 AU, Neptune 30x 1.561 1.698 1.762 1.772 

2 AU, Jupiter lx 1.367 1.329 1.305 1.281 

2 AU, Jupiter 3x 1.366 1.328 1.300 1.267 

2 AU, Neptune lOx 1.397 1.356 1.313 1.251 

2 AU, Neptune 30x 1.400 1.357 1.307 1.243 

5 AU, Jupiter lx 1.316 1.286 1.271 1.261 

5 AU, Jupiter 3x 1.297 1.270 1.257 1.248 

5 AU, Neptune lOx 1.523 1.529 1.463 1.335 

5 AU, Neptune 30x 1.523 1.530 1.462 1.349 

10 AU, Jupiter lx 1.372 1.312 1.277 1.253 

10 AU, Jupiter 3x 1.407 1.327 1.275 1.240 

10 AU, Neptune lOx 1.544 1.527 1.436 1.275 

10 AU, Neptune 30x 1.543 1.527 1.433 1.279 
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Table 6: Calculated Bond albedos (Ab) of the exoplanet model Jupiter and Neptunes in this work, 
integrating from 0.35 to 2.5 /im over geometric albedo spectra and phase integrals from the albedo 
spectral model described in ^3.2.51 and Appendix [XI Contributions from longer wavelengths than 
2.5 fim. to the Bond albedo are smaller than 5%. 



Model 


A B 


0.8 AU, Jupiter lx 


0.408 


0.8 AU, Jupiter 3x 


0.333 


0.8 AU, Neptune 10 x 


0.297 


0.8 AU, Neptune 30 x 


0.229 


2 AU, Jupiter lx 


0.914 


2 AU, Jupiter 3x 


0.899 


2 AU, Neptune 10 x 


0.787 


2 AU, Neptune 30 x 


0.685 


5 AU, Jupiter lx 


0.675 


5 AU, Jupiter 3x 


0.592 


5 AU, Neptune 10 x 


0.407 


5 AU, Neptune 30 x 


0.366 


10 AU, Jupiter 1 x 


0.489 


10 AU, Jupiter 3x 


0.348 


10 AU, Neptune 10 x 


0.378 


10 AU, Neptune 30 x 


0.344 



